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Abstract 

The idea of Media-based Modulation (MBM), introduced in 00 is based on embedding information in 
the variations of the transmission media (channel states). This is in contrast to legacy wireless systems where data 
is embedded in a Radio-Frequency (RF) source prior to the transmit antenna. MBM offers several advantages vs. 
legacy systems, including “additivity of information over multiple receive antennas”, and “inherent diversity 
over a static fading channel”. MBM is particularly suitable for transmitting high data rates using a single 
transmit and multiple receive antennas (Single Input-Multiple Output Media-Based Modulation, or SIMO- 
MBM). However, complexity issues limit the amount of data that can be embedded in channel states using a 
single transmit unit. To address this shortcoming, the current article introduces the idea of Layered Multiple 
Input-Multiple Output Media-Based Modulation (LMIMO-MBM). LMIMO-MBM enables forming a high-rate 
constellation as superposition of constituent vectors due to separate transmit units. Relying on such a layered 
structure, LMIMO-MBM can significantly reduce both hardware and algorithmic complexities, as well as the 
training overhead, vs. SIMO-MBM. Exploiting the proposed layered constellation structure, a fast iterative 
algorithm is proposed for signal detection, and a practical (small size and low complexity) RF configuration 
is presented for embedding information in channel states. Simulation results show excellent performance in 
terms of Symbol Error Rate (SER) vs. Signal-to-Noise Ratio (SNR). For example, a 4 x 16 LMIMO-MBM 
is capable of transmitting 32 bits of information per (complex) channel-use, with SER cx 10 -5 at E^/Nq cx 
— 3.5dB (or SER cx 10 -4 at E^/Nq = —4.5dB). This performance is achieved using a single transmission (no 
extension in time/frequency), and without adding any redundancy for Forward-Error-Correction (FEC). This 
means, in addition to its excellent SER vs. energy/rate performance, MBM relaxes the need for complex FEC 
structures used in legacy wireless systems, and thereby minimizes the transmission delay. Application of FEC 
can further improve the performance. For example, applying Reed-Solomon codes enables transmitting 30 bits 
of information per (complex) channel-use with a Frame Error Rate (FER) ~ 10“ 5 at E^/Nq ~ — 6dB. It 
is shown that, under a set of mild conditions, by applying FEC with error correction capability t , the slope 
of the error rate vs. SNR (with hard decision decoding) will asymptotically increase by a factor of t + 1. 
Overall, LMIMO-MBM provides a promising alternative to MIMO and Massive MIMO for the realization of 
5G wireless networks. 


I. Introduction 

Shannon capacity results indicate that the transmission rate can increase as a linear function of the 
available spectrum, multiplied by a logarithmic function of the transmit energy. Wireless communica¬ 
tions relies on two key attributes, which are also traditionally considered at its inherent bottlenecks. 
First, the spectrum is shared, resulting in mutual interference among wireless links operating over the 
same spectrum. Second, transmission channel includes a multitude of propagation paths, resulting in 
multi-path fading. Multi-path fading in turn can result in deep fades when signals received through 
different transmission paths add destructively. In many scenarios of practical interest, the transmission 
paths change only slowly with time (slow fading), potentially resulting in a long lasting degradation 
of the received Signal-to-Noise-Ratio (SNR), referred to as deep fades. 


Multiple-Input Multiple-Output (MIMO) antenna systems embrace the above two attributes towards 
improving the spectrum/power efficiency 000. as well as in dealing with deep fades In 
particular, the fact that signals transmitted from different transmit antennas will be superimposed at 
each receive antenna results in an input-output relationship captured in the form a matrix multiplication. 
The fact that such a channel matrix is, with high probability, non-singular is shown to (asymptotically) 
result in a linear scaling of rate with min (At,Ar), where At and Ar are the number of transmit 
antennas and the number of receive antennas, respectively. This phenomenon mimics the linear scaling 
of rate vs. spectrum by just using more transmit/receive antennas. It also enables combating deep 
fades by introducing redundancy among data to be transmitted/received through separate propagation 
paths. The improvement in spectrum efficiency, for large values of SNR, is captured in the so-called 
Multiplexing Gain, where MG= thvh{At , Ar) . MG appears as a scale factor, times the log (SNR), in 
the expression of rate as a function of SNR for SNR —> oo. It is also known that, to combat slow 
fading, the degrees of freedom offered by MIMO systems can be used to create diversity. Diversity 
order, which in essence captures the number of independent paths (independent gains) involved in the 
end-to-end transmission, can be increased, but only at the cost of a reduction in MG [7], 

Although MIMO systems provide an elegant way to tailor wireless communications to adopt to 
its two fundamental attributes/bottlenecks mentioned earlier, three issues limit their achievable rate 
vs. energy. First, the problem of deep fades still exists and can be only (partially) alleviated at the 
cost of a reduction in the achievable rate (MG) 0 Second, MG increases only with the smaller 
of the number of transmit and receive antennas. Third, the MIMO channel matrix is typically non- 
orthogonal, reducing the achievable rate as compared to an orthogonal channel matrix with a similar 
dimensionality. Media-based Modulation (MBM) deals with these three issues. The core idea is based 
on randomizing the wireless channel through perturbing the propagation environment in the vicinity 
of transmit antenna(s), which in turn will change the overall transmission path. This can be viewed as 
creating a multitude of channel states, each corresponding to a different transmission path, where the 
transmitter can select any of the channel states in each transmission. The information to be transmitted 
is used by the transmitter as an index to select a particular channel state. 

As mentioned above, MBM can be interpreted as creating a channel with a finite number of, say M, 
states. MBM transmitter selects one of these M states in each transmission, and thereby, can embed 
log 2 (M) bits of information in this selection. This can be realized if each channel state corresponds 
to a unique and selectable transmission path. Assuming the transmission paths corresponding to 
different channel states are independent of each other, for a static multi-path channel, the states 
correspond to independent realizations of the underlying channel gains, which are known to follow 
a complex Gaussian distribution. Later parts will discuss practical methods to independently change 
the transmission path for each channel state. Transmitter will then send a signal through the selected 
channel state, in order to convey the particular selection to the receiver. The process is such that the 
corresponding received signal acts as a unique signature for the selected channel state. This means, in 
the absence of noise, the received signal would uniquely determine the channel state. In the case of 
multi-path fading, the net effect is that each channel state results in a point over the receiver dimensions, 
and such points have independent, identically distributed (i.i.d.) Gaussian components. A key point 
is that, if K receive antennas are used, the number of complex receive dimensions will increase by 
K, while satisfying the i.i.d. property for the distribution of points over all the resulting coordinates. 
This means, constellation points will span the entire space formed over receiver K spatial dimensions. 
This phenomenon occurs even if a single antenna is used for transmission. This is in contrast to the 
traditional MIMO where the effective dimensionality is governed by the minimum of the number of 
transmit and receive antennas. The task of decoding at the receiver concerns mapping the signtuare 


signal, received in noise, to the corresponding channel state. Note that in this setup the primary role 
of the transmitted signal is to convey the signature of the channel state to the receiver, while in legacy 
transmission systems (called Source-Based Modulation, SBM, hereafter), the role of the transmitted 
signal is to convey information by selecting a particular candidate from a set of alphabets with a 
cardinality equal to the number of messages to be transmitted. It is well known that the rate associated 
with such an alphabet can be increased, but at the cost of an exponential increase in the transmit energy. 
MBM results in a new landscape for such an energy vs. rate trade-off. In MBM, the cardinality of the 
set of input signals is typically one, although one can embed further information by selecting the input 
signal from a set with a larger cardinality. This corresponds to combining MBM with legacy SBM. 

Example 1: Consider a wireless channel with two states si and s 2 . In one of the states, the channel 
gain is equal to 0.5, and in the other state, it is equal to 1.5. Transmitter can select either of the 
two states in each transmission, but is not aware of the correspondence of the two gain values to 
the two states (does now know which state corresponds to which gain value). Let us assume we are 
interested in transmitting two bits per channel use. In one scenario, reminiscent of SBM, the transmitter 
selects one of the two states for all its transmissions and use it with a Pulse-Amplitude-Modulation 
(PAM) constellation of size 4, composed of points {-3,-1,1,3}. With probability 1/2, the selected 
channel state corresponds to the one with the lower gain, resulting in the received constellation being 
composed of points C± = (—3/2, —1/2,1/2,3/2} with a d min = 1. This means one can guarantee a 
wrost case dmin equal to 1 using an average transmit energy of (9 + 1 + 1 + 9)/4 = 5. In a second 
scenario, reminiscent of MBM, for each transmission, transmitter uses one bit of information to select 
the channel state, and transmits a Binary-Phase-Shift-Keying (BPSK) modulation with points {—1,1} 
through the selected state. It easily follows that the receiver will observe a 4-PAM constellation with 
points C 2 = {—3/2,—1/2,1/2, 3/2}, again resulting in d m - m = 1, but this time at the cost of using one 
unit of energy vs. the 5 units used in the first scenario. The underlying assumption is that, although 
transmitter is oblivious to the actual realization of the channel gain, receiver is aware of the structure 
of the constellation and its labeling. This information can be easily conveyed to the receiver through 
an initial training phase, in which the transmitter selects the two channel states in an order that is 
prearranged with the receiver. For example, it is agreed to use s x in the first training period, and s 2 in 
the second one, or vice versa. It is also agreed that the state used in the first training symbol represents 
data bit zero and the second one represents data bit one, or vice versa. 

Although the above example is engineered to construct a constellation with regularly spaced points, 
in practice, the gains corresponding to different channel states correspond to realizations of a Gaussian 
random variable (multi-path fading in a rich scattering environment). As a result, the constellation 
formed at the receiver side will be composed of points with coordinates realized as samples of a 
Gaussian distribution. This is in accordance with capacity achieving random codebooks on AWGN 
channels. Although in practice traditional modulation schemes rely on uniformly placed points, sim¬ 
ulation results show that such a Gaussian random constellation typically performs well. A second 
example will follow to clarify this point. Indeed, reference [j2j shows that, as the constellation size 
increases, such a Gaussian constellation becomes capacity achieving. 

The above example shows another notable property. If the channel state is changed randomly from 
transmission to transmission, the equivalent channel can become Ergodic in time. In other words, the 
better and the worse channel states collectively contribute to constructing a single constellation, which 
will be in effect in all transmissions. This case corresponds to constellation C 2 in this example. In 
other words, low channel gains result in constellation points closer to the origin, while high channel 
gains result in constellation points further away from the origin. As a result, constellation points are 


spread across the available transmission space. This feature is analogues to traditional modulation 
constructions, such as PAM/QAM, which include constellation points over different energy shells 
distributed in the constellation space. This includes points near the origin (of low energy values) and 
far away from the origin (of high energy values), all contributing to the constellation cardinality and 
spacing. As a consequence, the inherent bottleneck occurring in the case of SBM, corresponding to 
being stuck with the worse channel state in all transmissions is avoided. This phenomenon is further 
discussed in later parts of the current article and is interpreted as an inherent (built-in) diversity effect. 
Reference [2] shows that, this phenomenon will asymptotically (for large constellation sizes) convert a 
static fading channel into an (Ergodic) AWGN channel. In other words, the harmful effect of multi-path 
fading, resulting in outage behaviours in deep fades, is (asymptotically) alleviated. The important point 
is that in MBM, this so-called built-in diversity is inherently (automatically) realized at no cost, in 
particular without wasting rate or energy. This is unlike traditional MIMO systems where an increase 
in the diversity order is unavoidably accompanied by a reduction is rate (MG). ■ 

Example 2: Consider a wireless system with a single transmit and a single receive antenna. The 
channel has 256 states, each resulting in a constellation point with independent, identically distributed 
(i.i.d.) Gaussian components over the single (complex) receive dimension. Each realization of such 256 
points has an achievable rate (assuming AWGN and equal probability for the points). The achievable 
rate is a random variable depending on the particular realization of points. Figure [T] shows such 
achievable rate values for 10000 realizations, as well as the achievable rate of a 256 Quadrature 
Amplitude Modulation (QAM) constellation (again, assuming AWGN and equal probability for the 
points). It is observed that the distribution of achievable rate values deviates only slightly from that of 
the 256QAM. This means, with a slight margin for energy (in the order of ldB for the SNR point shown 
in details in Figure [TJ, the outage probability for the random constellation can become negligible. This 
significantly outperforms the outage behaviour of a 256QAM over a static Raleigh fading channel. 
Indeed, as the constellation size increases, the capacity tends to the capacity of an (Ergodic) AWGN 
channel and the underlying outage probability will tend to zero [2]. This means, unlike using a uniform 
constellation with points occurring with equal probabilities, the random Gaussian constellation will also 
realize the shaping gain. ■ 

A. Literature Survery 

Reference [1] shows that embedding part or all of the information in the (intentional) variations of 
the transmission media (channel states) can offer significant performance gains vs. traditional Single- 
Input Single-Output (SISO), Single-Input Multiple-Output (SIMO) and Multiple-Input Multiple-Output 
(MIMO) systems. This method, coined in [1] as Media-Based Modulation (MBM), is in contrast with 
traditional wireless systems where data is embedded in the variations of an RF source (prior to the 
transmit antenna) to propagate via fixed propagation paths (media) to the destination. In particular, 
using capacity arguments, reference 0 shows that by using a single transmit antenna and a single or 
multiple receive antennas; MBM can significantly outperform SBM. 

Following 0. reference [2] proves that, alx I\ MBM over a static multi-path channel asymptotically 
achieves the capacity of K (complex) AWGN channels, where for each unit of energy over the single 
transmit antenna, the effective energy for each of the K AWGN channels is the statistical average of 
channel fading. In addition, the rate of convergence is computed. It is shown that significant gains 
can be realized even in a SISO-MBM setup. An example for the practical realization of the system 
using RF mirrors, accompanied with realistic RF and ray tracing simulations, are presented. Issues of 
equalization and selection gain are also briefly discussed. 
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Fig. 1: Capacity of random vs. regular 256 points constellation. 


The idea of embedding information in the state of a communications channel is not new. MachZehn- 
der modulators, widely used for signalling over fiber, modify the light beam after leaving the laser. 
However, due to the lack of multi-path in single mode fibers, the advantages due to SIMO-MBM and 
MIMO-MBM, realized in the context of wireless, do not apply. 

In distant relationship to MBM, there have been some recent works on embedding data in antenna 
beam-patterns |(8j-[10] or antenna selection [ 1_1 ]-[ 15]. Note that unlike MBM, none of these works can 
realize advantages due to embedding information in the channel state. Most notably, these advantages, 
reported for the first time in [1] |2], include “additivity of information over multiple receive antennas” 
and “inherent diversity without sacrificing transmission rate”. 

In 0 data is embedded in two orthogonal antenna beam-patterns, which can transmit a binary sig¬ 
nal set. Although use of orthogonal basis is common in various formulations involving communications 
systems, it usually does not bring any benefits on its own, it just simplifies problem formulation and 
signal detection by keeping the noise projections uncorrelated. This means there are no clear advantages 
in designing the RF front-end to support orthogonal patterns as used in [8] j9}. The motivation in j8] 
m is to reduce the number of transmit chains and no other benefits are discussed. 

Bains |T0j discusses using parasitic elements for data modulation, and shows limited energy saving, 
which again is due to the effect of classical RF beam-forming. 

Spatial Modulation (SM) [11]-[15] uses multiple transmit antennas with a single RF chain, where 
a single transmit antenna is selected according to the input data (the rest of the data modulates the 
signal transmitted through the selected antenna). SM is in essence a diagonal space-time code, where 















the trade-off between diversity and multiplexing gain has been in favour of the latter. A shortcoming 
of SM is that the rate due to the spatial portion increases with log 2 of the number of antennas, while in 
MBM, it increases linearly with the number of RF mirrors (on-off RF mirrors are introduced in [j2] as 
means of embedding binary data in the channel state). In SM, antennas should be sufficiently separated 
to have independent fading, while in MBM, RF mirrors are placed side by side. The switches used in 
SM are high power, which means expensive/slow, or each antenna needs a separate Power Amplifier 
(PA) with switches placed before PAs. The switches used for RF mirrors in MBM are cheap, low 
power and fast. 

The use of tuneable parasitic elements external to the antenna(s) for the purpose of RF beam-forming 
is well established. However, the objective in traditional RF beam-forming is “to focus/steer the energy 
beam, which does not realize the advantages of MBM (where data is modulated by tuning external 
parasitic elements). 

In continuation to SM [11|-[[T5|, Space-Shift Keying (SSK) [ 16]-[ 18] and Generalized SSK (GSSK) 
[19] have been studied for low-complexity implementation of MIMO systems. Again, the key motiva¬ 
tion behind the application SM/SSK/GSSK in [ 1 1]-|[T9| is the use of a single RF chain, and accordingly, 
one antenna remains active during data transmission. Other advantages are mentioned as avoiding 
inter-antenna synchronization and removing inter-channel interference JT8| ]. In addition to complexity 
considerations, it is shown that these modulation schemes may offer better error performance as 
compared to conventional MIMO techniques [ [T6| . 

The advantages of MBM, which are discussed in details in [iJQ, are briefly explained next. 


B. Advantages of MBM 

MBM offers an inherent diversity in dealing with static fading. Since constellation points in MBM 
correspond to different channel realizations, the spacing among points is formed using different channel 
gains. Hence, both good (high gain) and bad (low gain) channel conditions contribute towards forming 
the required spacing among constellation points. This feature of MBM removes the bottleneck of deep 
fades in static fading channels without compromising the rate. This is unlike traditional MIMO setups, 
where an increase in the diversity order is accompanied by a decrease in rate (multiplexing gain). It 
was shown that significant gains can be realized even in a SISO-MBM setup due to inherent diversity 
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In 1 x K SIMO-MBM, due to the independence of channel gains to different receive antennas, the 
received signal constellation spans all the Q = 2 x K receive dimensions. Therefore, MBM benefits 
from larger spacing among constellation points in case of using multiple receive antennas. This is 
analogous to additivity of information over multiple receive antennas. This is in contrast to SIMO- 
SBM where the received signal spans a single complex dimension, and consequently, only SNR gains 
can be achieved through techniques such as maximum ratio combining (without increasing the effective 
signal space dimensionality). In addition, a desirable feature of MBM is that the constellation points 
are spread across the signal space following a Gaussian distribution (due to Raleigh fading), which 
is in agreement with Gaussian random coding requirements. Relying on this feature, our earlier work 
reported in (2) shows that a 1 x K MBM over a static multi-path fading channel asymptotically achieves 
the capacity of K parallel (complex) AWGN channels, where for each unit of energy over the single 
transmit antenna, the effective energy for each of the K AWGN channels is the statistical average of 
the channel fading. 

Another benefit of SIMO-MBM is the “K times energy harvesting” which means, assuming K 
receive antennas and a fading with statistical average gain of one, the average received signal energy 





is K times the transmit energy. Legacy MIMO systems enjoy a similar property, however, in the case 
of MIMO systems, the channel matrix is typically non-orthogonal. This results in correlation among 
noise components along different receive dimensions if channel inversion is used at the receiver. Eigen 
beam-forming can be used to diagonalize the channel matrix, but the underlying issue will surface in 
another equivalent form; it results in different channel gains along different eigen-dimensions. Both of 
the above phenomena (caused by the channel matrix being non-orthogonal) will reduce the capacity 
vs. a set of parallel channels with independent noise components and equal SNR. In contrast, in MBM, 
noise components are added at different receive antennas and are independent of each other and have 
equal variances. This results in equal SNR for all the underlying parallel channels, and increases the 
capacity as compared to legacy MIMO. To further explain this phenomenon, note that in applying 
water filling to achieve the capacity of a K x K MIMO channel, for low energy values, only one 
eigenvector receives all the energy. Although this is the eigenvector with the largest eigenvalue, it still 
means a channel with a single dimension (rather than K dimensions which will be achieved only if 
the energy level is high enough such that all eigenvectors are occupied). Readers are referred to 0 
for further discussions. 

Some complexity issues arise when only a single RF transmit unit is used to realize the advantages 
of MBM at high data rates. Next section discusses the underlying practical issues, and present methods 
to address them. 

C. Limitations of SIMO-MBM in Transmitting High Date Rates 

Using a single transmit unit to embed all the information in channel states will limit the amount 
of information that can be practically transmitted per channel-use. For example, let us assume we are 
interested to transmit 32 bits of information per channel-use. The complexity in using a single RF 
transmit unit to encode all 32 bits may be excessive. The reasons are: 

C.l. It is practically difficult to use 32 RF mirrors in a single RF transmit unit (on-off RF mirrors 
are introduced in [2] as means of embedding binary data in channel states). 

C.2. Training requires transmitting 2 32 test signals, which is resource intensive and is vulnerable to 
channel time variations. 

C.3. Detection requires searching (minimum distance decoding) among 2 32 signal points, resulting 
in excessive algorithmic and storage complexities. 

C. 4. To deal with channel time variations, it is of interest to track the changes in the position of the 

constellation points in order to increase the minimum time interval between successive training 
phases. It is difficult to track 2 32 constellation points. 

D. Proposed Solution 

To address the above issues, this article proposes a new method to perturb the RF channel, resulting 
in a layered constellation structure. Assume several separate transmitter units, each generating a set 
of received vectors (called “constituent vectors” hereafter), operate at the same time. As a result, the 
received vector will be the sum of the constituent vectors corresponding to different transmit units. 
As the constituent vectors are random and independent from each other, the cardinality of the set of 
received vectors will be (with high probability) equal to the product of the cardinalities of the set of 
constituent vectors corresponding to different transmit units. As a result, using R n RF mirrors at N 
transmit units (each unit has a single radiating element) results in 2 NxRn received vectors, capable 
of transmitting R = N x R n bits of information per channel-use. Transmit units are arranged such 
that there is a negligible coupling among them. As a result, constellation points will be formed as the 


sum (superposition) of constituent vectors due to each transmit unit. To emphasize this superposition 
property, which forms the basis behind the complexity reduction, the proposed approach is called 
Layered MIMO-MBM, or LMIMO-MBM, hereafter. Benefits of this setup in reducing complexity are 
explained next. 

D.l. Number of RF mirrors used at individual transmit units is reduced by a factor of N. 

D.2. Detection is performed using an iterative search algorithm. At each step, the proposed algorithm 
searches for the constituent vector contributed by a given transmit unit (initializing the search 
by zero vectors) and continues iteratively. To improve the search result, one can start from 
multiple initial points (for example, corresponding to different permutations of transmit units) 
and at the end select the best candidate. 

D.3. Training is simplified as it is composed of N separate training tasks, each over a smaller set of 
alphabet size 2 Rn , as compared to training over 2 NxRn elements. 

D.4. Tracking is simplified as it is composed of N separate tracking tasks, each over a smaller set 
of alphabet size 2 Rn elements, as compared to tracking 2 NxRn elements. 


For example, to send 32 bits of data per channel-use, one can use 4 transmit units each modulating 
8 bits, which means only 8 on-off RF mirrors are required in each transmit unit. Training/tracking is 
composed of 4 separate tasks, each involving a small alphabet size of 2 8 = 256 elements. 

A disadvantage of the LMIMO-MBM method proposed in the current article is that the received 
constellation vectors do not correspond to independent Gaussian vectors any longer as the constellation 
vectors are summation of a smaller number of independent constituent (Gaussian) vectors. This is 
in contrast with the requirement of Gaussian random coding over AWGN channels. However, the 
corresponding degradation in SNR performance is negligible. For example, Figure [2] demonstrates an 
example for this gap when minimum distance decoding is performed using exhaustive search. Since 
exhaustive search is not feasible for larger constellations, Figure [5] shows the performance of the 
proposed sub-optimal iterative decoder (to be explained later) for a constellation size of 2 32 . 

In the following, first the system model for SIMO-MBM and LMIMO-MBM are reviewed in section 
[II| The benefits of using multiple transmit units in LMIMO-MBM are studied in more details in 
section [nl| Subsequently, in section [TV} a sub-optimum iterative algorithm that enables fast detection 
in LMIMO-MBM setup is proposed. In section |V} the LMIMO-MBM performance curve and the 
simulation setup are discussed. Section [VI] discusses the application of FEC to LMIMO-MBM, includ¬ 
ing some simulation results. It is shown that, under a set of mild conditions, by applying FEC with 
error correction capability t, the slope of the error rate vs. SNR (with hard decision decoding) will 
asymptotically increase by a factor of t + 1. The issue of MBM potential bandwidth increase, due its 
time varying nature, is addressed in section |VII[ Finally, a practical (small size and low complexity) 


RF configuration for embedding information in channel states is presented in section VIII 


II. System Model 

A. SIMO-MBM 

In a 1 xK SIMO-MBM system (see Figure [3]), where there are M messages to be sent, each message 
selects a channel realization h(m) with complex components ///, (rn), k = 1..... A' and rn = 1,..., M. 
E| h k {m) | 2 = 1 where E denotes statistical averaging. Additional bits of information, corresponding 
to SBM message s, can be transmitted by directly modulating the RF signal. As an example, this can 
be done using RF phase shifters to embed the SBM information bits in the phase of the RF signal 
(e.g., selecting O' or 180' phase shifts to generate S(s) for sending one additional bit). Assuming linear 
modulation for the SBM part, the SBM message s can be viewed as a complex number multiplying 









Fig. 2: Performance of the 2x8 LMIMO-MBM (2 RF mirrors are used in each transmit unit) vs. 1x8 
SIMO-MBM (4 RF mirrors are used in a single transmit unit), R = 16 bits per channel-use for a single 
transmission in time (without any Forward-Error-Correction). Note that at the cost of small degradation 
in performance and increasing the number of transmit units from one to two, the 2x8 LMIMO-MBM 
offers reduced complexity in each RF structure and reduction in algorithmic complexity for training, 
detection and tracking. 


the RF carrier prior to the antenna, which will also appear as a scale factor in the received signal. In 
other words, the received signal (prior to AWGN) will be equal to sh(m). To simplify presentation, 
the dependency on SBM message s is ignored in what follows. AWGN z has independent identically 
distributed (i.i.d.) components z k , k = 1 ,..., K, E| z k | 2 = N 0 . In the absence of Forward-Error- 
Correction (FEC), the total transmission rate is equal to R = log 2 M bits per channel-use. The receiver 
is trained to know h(m). This is done by sending a separate training signal for each possible channel 
realization (state) to measure received constellations h k (m), 'ini. ik. The order of scanning through 
channel states in the training phase is prearranged such that the assignment of channel states to data 
indices is implicit. For example, the channel state indexed by data m is selected during the mth 
training symbol (to covey h(m) to the receiver). For a Rayleigh fading channel, h k (m),k = 1,..., K, 
m = 1, are i.i.d. complex Gaussian. 

B. Layered MIMO-MBM (LMIMO-MBM) 

In a N x K MIMO-MBM system (see Figure [4]), the M messages to be sent are distributed 
among N transmit units, using vector m with components mi,.... rn^. Subsequently, each transmit 















































rn : MBM message, m£|l,...,2* m | 
S : SBM message, S £ |l,..., 2 Rs | 

R = R m +R s 



S(s) , specifies the combination 
of Baseband and RF modulation 


Fig. 3: SIMO-MBM block diagram (assuming linear modulation for SMB data). 


unit selects its own channel realization vector h n (m n ), n = according to m n . These channel 

realization vectors, called constituent vectors, add to form the received constellation points c(m). Again, 
E| h'^(m n ) | 2 = 1, Vm, Vn and Vfc. Similar to SIMO-MBM, each transmit unit can send optional SBM 
information bits. To attain this, transmit unit n,n = 1. .... N, modulates its own RF signal based on 
the nth component of the information vector s, to generate RF modulated signal S n (s n ), n = 1,..., N. 
Due to symmetry, distributing rate and power equally among individual transmit units results in the 
highest end to end mutual information. Throughout the paper, we assume power/rate is equal among 
transmit units. Consequently, the sets of LMIMO-MBM constituent vectors, corresponding to different 
transmit units, are of equal cardinalities, namely \f~M = 2 Rrn ^ N = 2 Rn . Assuming: (1) unit power for 
each transmit unit, and (2) linear modulation for the SBM part, the received signal will be equal to: 

N 

c(m) = J^s n h n (m n ), (1) 

n=l 

where c(m) is the received constellation point with K complex dimensions. This means constellation 
points are formed as a weighted superposition of constituent vectors due to each transmit unit. For a 
Rayleigh fading channel, elements of channel realization vectors h^(m n ) are i.i.d. complex Gaussian. 
However, the received constellation vectors are no longer independent. Again, to simplify presentation, 




Fig. 4: LMIMO-MBM block diagram (assuming linear modulation for SMB data). 


the dependency on SBM message is ignored in what follows. As a result, [T] is simplified to: 

N 

c(m) = J^h n (m n ). (2) 

71=1 


III. Benefits of Using Multiple Transmit Units in LMIMO-MBM 

Using multiple transmit units simplifies receiver training and signal detection, specifically in trans¬ 
mitting large amounts of information per channel-use. Although MBM transmitter is oblivious to the 
details of the constellation structure, receiver training is required to convey the constellation structure 
to the receiver side. Training complexity grows linearly with the size of the constellation, making it an 
expensive task for constellations of large cardinalities. When using multiple transmit units in LMIMO- 
MBM, constellation set is formed by superposition of smaller sets of size vM. Therefore, one can 
measure V«, Mk and Vm, by sending pilot signals (using one transmit unit at a time) and then 

constructing the entire constellation at the receive side by using superposition property (captured in 
equation ([!])). This means complexity of the receiver training will reduce from O(M) to 0(N X y/M). 
Furthermore, there would be no need for huge memory to save all the constellation points c, since 
on demand computation can take place, provided that constituent vectors h n are known at the receive 
side. 






Another benefit is the reduction in the complexity of signal detection by exploiting the structure 
of the constellation due to its layered construction. In this article, this benefit is realized using a 
(greedy/sub-optimum) iterative search algorithm, which searches among \/M constituent vectors at 
each iteration. 

In the proposed LMIMO-MBM, the dependency among constellation points results in a small 
decrease in the end-to-end mutual information, as compared to using random constellations, i.e., a 
constellation with i.i.d. Gaussian distribution for the components of the constellation points (resembling 
random Gaussian code-books). However, the performance degradation remains negligible (for example, 
see Figure [2]). 


IV. Detection Algorithm 

In the presence of AWGN, for a given received signal r = c + z, maximum likelihood detection 
for constellation point c corresponds to c, which is at minimum Euclidean distance to r. A naive, 
but optimal solution for finding c is exhaustive search among all constellation points. This so-called 
“linear search” is guaranteed to find the nearest point with a complexity of O(M). However, this 
operation can be very expensive for large signal sets. The problem of finding the nearest member (in 
our case, in Euclidean distance sense) of a set of predefined points to a query point, also known as 
Nearest Neighbor Search (NNS), has many applications, including in pattern recognition, statistical 
classification, computer vision, etc. Therefore, the NNS problem has been widely studied and many 
algorithms have been suggested in the literature. However, none of the existing classes of algorithms 
appeared to be suitable for LMIMO-MBM setup due to the reasons explained next. 


A. A Review of Potentially Suitable Search Algorithms 


A general approach for reducing the complexity of a full search in NNS problem is based on 
preprocessing data points (in our case, constellation points) and storing them in a suitable data structure, 
such as k-d trees (20) . This provides a near optimal solution for cases of low dimensionality. However, 
the complexity of most these algorithms grow exponentially with the number of dimensions, which 
is the case when multiple receive antennas are used in MBM. This results in a search complexity 
comparable to that of the naive “linear search”. Furthermore, due to the nature of propagation media, 
channel realizations, and hence the constellation points, change in time. Therefore, a detection algorithm 
which relies on extensive preprocessing of the constellation points will not be favorable. Another 
approach to reduce the complexity is based on Triangle Inequality Elimination (TIE) tf2l| - (24). In this 
class of methods, reference points are used and their distances to actual data points are precomputed 
and stored. In the search phase, the distances between the query and the reference points are calculated. 
Then, using the precomputed distances and a simple comparison, a large number of data points are 
eliminated. A drawback of this technique is its large memory requirement (larger than the constellation 
size itself) to store the precomputed distances, as well as the computational overhead due to the required 
preprocessing. 

In the following, a fast greedy iterative algorithm, enabled by layered structure of the proposed 
MIMO-MBM, is used for sub-optimal minimum distance signal detection. Advantages of the proposed 
algorithm include: 


• The complexity grows as 0( \/M) compared to 0{M) in linear search. 

• The complexity grows linearly with the number of dimensions in contrast to the alternative 
search mechanisms known in the literature of NNS. 



• The algorithm does not require any storage space more than what is required to store channel 
realizations (constituent vectors) for the individual transmit units, which is of order 0(Nx y/~M). 

• There is no need for preprocessing, in contrast to alternative search mechanisms known in the 
literature of NNS. 


B. MIMO-MBM Greedy/Iterative Search Algorithm 

Iterations: Algorithm [T] is the pseudocode for the proposed sub-optimum minimum distance detec¬ 
tion in LMIMO-MBM. Each iteration is composed of N greedy steps. At each step, the algorithm 
attempts to (successively) find the constituent vector corresponding to a single transmit unit. At step 
i, a potential solution c (l) , for the minimum distance vector to the query signal r is given as the 
current best solution. Suppose is formed by superposition of constituent vectors (h 1 , ...,h iV ). The 
algorithm compares the constituent vector corresponding to transmit unit i, which means h\ with all 
other vectors in the set W = { IT( 1),..., IT ( 2 Rn )} of constituent vectors corresponding to transmit unit 
i. Subsequently, IT is replaced with the vector IT in W that results in a constellation point c (l+1} , 
which (among all choices available at this step) is the closest to the signal r. Therefore, 


argmin| r — (c^ — h* + h) | 2 

(3) 

h ew 


c (i+1) = c (i) -tf + hL 

(4) 


In other words, at each step, a search over constituent vectors corresponding to a single transmit unit 
is performed (based on the current knowledge of signals of the other transmit units). Each iteration is 
composed of N such steps, and in total T iterations are performed before concluding the search. Upon 
completion of each iteration, c( N+1 '> is passed to the next iteration as the current best candidate, i.e., 
closest to the query signal r. As can be seen in ([3]), there is no need to store all constellation points, 
since they can be computed on the fly using constituent vectors. 

Initialization: At the first step of the first iteration, there is no prior knowledge of any of the 
constituent vectors. Therefore, c ( 1 ’ is initialized to a vector of all zero elements, which in fact is not 
a valid constellation point. At the end of the first iteration (i.e., after N steps), c ( " jV+1) will be a valid 
constellation point. 

Improving the Search: The algorithm in its simple form explained above can find the optimum 
solution in vast majority of cases. However, since we are dealing with very low error rates, even a 
small number of errors due to the algorithmic failure (e.g. iterations end-up in a loop) can govern 
the performance and result in error floor. Simulation results show that, given a received signal r, the 
decoded constellation point may be different when different orderings of transmit units for (successive) 
decoding are considered. For this reason, multiple runs of the above algorithm, corresponding to 
different permutations of transmit units, is used to improve the search. Each of these runs correspond 
to one of AH possible permutations of the sequence (1,..., N), corresponding to a different ordering 
of transmit units in successive decoding. 

The performance can be further improved by using multiple (i.e P > 1) candidate points in each 
run. This means at step i, given P best points as the current solution, the algorithm replaces the Ah 
constituent vector in each of these P points with constituent vectors in the set W, and updates the set 
of the best P points for the next step (or next iteration). Larger values for P increases the chance of 
finding the nearest point to r. 



Algorithm 1 Iterative detection algorithm 


Search for the nearest (Euclidean distance) constellation point to received signal r 

1 

function FiND(r) 


2 

C( JV +!) = 0 


3 

for j = 1,..., T do 

D> j is the loop index for iterations 

4 

c(!) = q(N+1) 

5 

for i = 1,..., N do 

o i is the loop index for steps 

6 

h* = argmin r — (c^ — h* + h) | 2 

hew 


7 

c (i+ i) = cW - IT + h* 


8 

end for 


9 

end for 


10 

return 

> nearest point 

11 

end function 



V. Simulation Setup 

For the simulation of the proposed LMIMO-MBM system, the constituent vectors corresponding to 
each transmit unit are generated with complex i.i.d. Gaussian random components of unit variance, 
and then their linear combination is used to form the constellation. The performance is averaged over 
many independent runs. E b , energy per bit, is defined as the sum of signal energies of transmit units 
divided by the total number of bits per transmission, and N 0 is the AWGN spectral density at individual 
receive antennas. Therefore, 


E b /N 0 = 


NE 

RNo 


Es 

RN 0 J 


(5) 


where E is the signal energy corresponding to each transmit unit, and E s is the total transmit energy 
per constellation point. 

Figure [5] shows the SER performance for 4 x 16 LMIMO-MBM, over a static Rayleigh fading 
channel with AWGN, transmitting 32 bits of raw data per channel use. Detection is concluded after 
two iterations (over 4 transmit units), using P = 128 candidate points and considering different 
permutations of transmit units. 

Note that the performance shown in Figure [2] and Figure [5] is over a single transmission, without 
using any FEC. It can be argued that in MBM, there is no need to use complex FEC structures, which 
are typical in legacy wireless systems. In the absence of FEC, signal detection can be performed with 
the lowest possible delay of a single symbol. These features of “lowest possible decoding delay” and 
“small error probability over a single time transmission” simplify the use of methods based on decision 
feedback (there will be no error propagation effects), for example for the purpose of equalization and/or 
tracking the constellation structure in time. However, the addition of FEC to MBM is straightforward 
and can be realized with simple code structures with hard decision decoding, such as BCH (Bose- 
Chaudhuri-Hocquenghem) codes or RS (Reed-Solomon) codes. This issue will be discussed next. 


VI. Application of Forward Error Correction Codes 

In the application of FEC to MBM, it is beneficial to consider coding schemes that operate on MBM 
symbols. Class of Group codes with alphabet size M would be a natural choice. Reference [27] proves 











Fig. 5: Performance of 4x16 MIMO-MBM for a single transmission (without any FEC), R = 32 
bits per channel-use, detection is performed using Algorithm [T] considering different permutations of 
transmit units with P = 128, T = 2. 


that in searching for good group codes, one can limit the search to those formed over elementary abelian 
groups. Maximum Distance Separable (MDS) codes, including Reed-Solomon codes, are a subclass 
of such group codes for which the minimum distance d min has the maximum possible value satisfying 
the Singleton bound. As a result, a natural choice in the application of FEC to MBM would be the 
class of Reed-Solomon codes with alphabet size M and block length M — 1. Reed-Solomon codes 
used in this article are obtained by puncturing such a larger code. Using such a Reed-Solomon code 
of block size L and dimension D with minimum distance d min results in an error correction capability 
of t = . 

Figure [6] shows the performance in the application of Reed-Solomon codes with symbols correspond¬ 
ing to MBM constellation points. Figure [7] shows the FER vs. SER performance of such Reed-Solomon 
codes. 

Appendix A shows that by applying such a coding scheme in conjunction with hard decision 
decoding, the slope of the error rate vs. SNR will asymptotically increase by a factor of t + 1. 
Consequently, application of FEC will realize an effect similar to “diversity over time” with a diversity 
order determined by t + 1. 
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Fig. 6: Performance of 4x16 MIMO-MBM, uncoded vs. Reed-Solomon code. 


Note that, unlike random-like codes such as Turbo-code and Low Density Parity (LDPC) codes 
which typically suffer from error floor, the slope of the error curve in coded MBM will not change as 
SNR increases. 

VII. Pulse Shaping, Equalization and Decision Feedback at the Receiver 

MBM inherently corresponds to a linear time varying system. Unlike linear time invariant systems, 
which maintain the spectrum occupancy of the signal, a linear time varying system can potentially 
increase the bandwidth. The issue of increase in band-width can be tackled by using the combination 
of the two methods explained next. 

Pulse Shaping: This corresponds to applying a pulse shaping filter at the transmitter in the time 
domain, such as raised cosine, to smoothen the transition between successive transmission pulses. 






Fig. 7: Frame Error Rate vs. Symbol Error Rate performance of 4x16 MIMO-MBM with Reed- 
Solomon codes. 


Silent Transmissions: This corresponds to inserting gaps of zero transmission (transmitter remains 
silent) between successive transmission pulses (see Figure [8]), and sampling at the receiver side 
corresponding to both on and off transmission intervals. Due to channel impulse response, the receiver 
will receive values even during the periods that the transmitter has remained silent. Receiver relies on all 
received values for signal detection. Note that impulse responses, corresponding to different channel 
states in MBM, are independent from each other. As a result, a longer channel impulse response 
increases the effective dimensionality of the receive signal space. For example, using D = 2 time slots 
in conjunction with 8 receive antennas results in 16 receive dimensions. Similarly, D — 4 time slots 
in conjunction with 4 receive antennas, result in 16 receive dimensions. Note that, if the length of the 
channel impulse response is larger than D, then there will be propagation from each pulse into its 
subsequent time intervals beyond the window of length D. Since we are dealing with very low SER 
values, the effect of this propagation beyond the window of length D can be removed using decision 
feedback. This means, receiver can decide for the transmitted message (corresponding to a given 
constellation in time) and account for its propagation into the next set of values (by subtracting the 
effect of the propagated signal). This is performed prior to the decoding of the subsequent constellation 
symbols. Simulation results show that, even for D = 2 which is realizable in a channel with an impulse 
response of a length greater than or equal to two, using raised cosine pulse shaping, the out-of-band 
radiation is significantly below the levels allowed in current wireless standards. Note that using D = 2 
allows for a bandwidth expansion by a factor of two, which is quite larger than the spectral expansion 
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are use to select channel states 


Fig. 8: Pulse shaping in time. 




Fig. 9: Two transmit unit stacked vertically. Note that more units can be stacked using a similar 
structure. 


due to time varying nature of MBM. 

As a side benefit, increasing dimensionality of the received signal space through silent transmissions 
allows reducing the number of receive antennas proportionally. 
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Fig. 10: Structure of a single RF mirror, all diodes are turned on and off (switched) simultaneously. 


VIII. RF Implementation 

Figure [9] shows 2 RF transmit units stacked uniformly along the z direction. Vertical stacking of 
transmit units reduces the RF coupling among them. The gap between adjacent units is 100 mm. 
Simulation result (using High Frequency Structural Simulator, HFSS) shows that this spacing results 
in about —25 dB coupling at 5.8 GHz. Each unit is formed using 14 RF mirrors. Printed Circuit 
Boards (PCBs) have a thickness of 60 mil. The height and diameter of the units are 50.75 mm and 
116.92 mm, respectively. Each PCB consists of two main parts with a spacing of 0.82 mm: (i) a 
monolithic rectangular copper part of size 46.75 mm x 12.5 mm, and (ii) seven rectangular patches 
with a length of 7.25 mm and width of 11.20 mm.The dimensions of patches are optimized at 5.8 GHz 
using HFSS to achieve high transmission efficiency (S tl < —15 dB) and low coupling (S tJ < —25 dB). 
Each two adjacent patches can be connected or disconnected using a diode (as an RF switch). In each 
RF mirror, all diodes are turned on or off (switched) simultaneously. The gap between two adjacent 
patches is 0.65 mm. A dipole antenna is placed at the center of each unit as the radiating element. 
Each RF mirror will be transparent to, or will reflect the incident wave, if its diodes are open or short, 
respectively. Each RF mirror in each transmit unit can be switched independently. Using 14 mirrors 
results in 2 14 channel states, where a subset of a smaller size, say 256, is used to transmit individual 
constituent vectors in the proposed LMIMO-MBM structure. 

To avoid using multiple transmit chains, a single transmit chain, i.e., single base-band, single RF 






















modulator, single Power Amplifier (PA), can be used as shown in Figure [TTJ The RF phase shifters 
can be used for sending additional (source-based modulated) data. For example, one can select 0’ or 
180’ phase shifts according to one additional bit of information (per transmit unit), or O’,90’, 180’,270’ 
phase shifts according to two additional bits of information (per transmit unit). The RF phase shifters 
can also facilitate training as will be explained next. 

In the training phase, to train for each transmit antenna, one possible option is to bypass (disable) the 
remaining antennas using RF switches (shown as optional components in Figure [TTJ). Another option 
is to train for each antenna when the other antennas are configured to transmit a default signal (e.g., 
transmit the constituent vector indexed by zero). In this case, the receiver will be first trained to learn 
the received vector corresponding to the transmission of the default constituent vector for each transmit 
unit. To train the receiver for the default constituent vectors, one can use the O’, 180’ phase shifts to 
realize ±1 towards forming of a Hadamard basis over transmit units. As a result, there is no need for 
RF switches in Figure [TTJ and trainings for default constituent vectors can be realized using phase 
shifters with at least two selectable phase values, 0’ and 180’, to multiply the signal from each transmit 
antenna by +1 or —1, respectively. In this case, the receiver will be able to invert the Hadamard basis 
and extract the received point corresponding to the default constituent vector for each transmit unit. 
Once the received points corresponding to default constituent vectors are known, the receiver will be 
trained by scanning through different constituent vectors corresponding to each transmit unit, while 
the other transmit units are set to send their default constituent vectors. Receiver will measure the 
signal received corresponding to the transmit unit being trained. Then, the receiver can compute the 
received constituent vectors due to any given transmit unit by accounting for the default constituent 
vectors corresponding to other units. 
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Fig. 11: Overall structure of a transmitter with multiple transmit antennas. 
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Appendix A 

Reference [27] defines the concept of “information set” as a subset of positions A within the coded 
block if, in the projectiorQ of the code over A, every possible 77-tuple of elements (77-fold Cartesian 
product alphabet) occurs exactly once. Reference [|27j proves that for MDS codes any subset of indices 
of size 77 is an information set. 

Let us assume there are M D code-words of length L forming a group code over an elementary 
abelian group of alphabet size M. Consider A < 77 positions within the FEC block. Let us project the 
code over these A positions. Given the alphabet size M, there can be at most M A different vectors 
of length A. Obviously, for A < 77, the set of vectors obtained through such a projection should have 
repeated elements. The underlying group property causes a group partitioning, and consequently, each 
combination is repeated an equal number of times, say V. It is straightforward to show that in general 
V > M D ~ A , while for good codes V = M D ~ A . 

To explain the above point further, let us focus on linear binary codes, i.e., M — 2, which are the 
most common form of elementary abelian group codes. It follows that the set of all code-words with 
0 in any given position form a sub-group, and the set of all code-words with 1 in that position form 
its coset. The number of code-words in the sub-group is equal to the number of code-words in the 
coset, unless all code-words belong to the sub-group, i.e., all code-words have a 0 in the position under 
consideration. A similar argument applies when instead of a single position, we consider a subset of 
positions V composed of A indices. The set of code-words with all zeros in those A positions form 
a sub-group, and there are up to 2 A — 1 cosets of this sub-group which are all of equal size. The 
number of cosets will be less than 2 A — 1, only if some of the combinations of bits do not occur in the 
projection of the code over 77. Considering all the 2 D binary code-words, this entails, in the projection 
of code over 77, either “all the possible 2 A combinations have equal frequency”, or “some combinations 
do not occur at all, which in turn means some combinations occur more often”. Good codes satisfy 
the former condition. In the language of information theory, this property means that any vector of 
size A, regardless of where the A positions are within the coded block, contains the same amount of 
information about the coded message, and this information has its maximum possible value of A bits 
for A < 77, and 77 bits for 77 < A < L. We refer to this property as the “Maximum Information 
Property” (MIP). In other words, all the positions within the coded block contribute equally to the 
information content of the code-words and their associated <7 min . If this property is not satisfied, the 
underlying code can be enhanced. This conclusion follows noting that some of the combinations in the 
underlying Cartesian product alphabet corresponding to indices in 77 do not occur in any code-word, 
and consequently, can be pruned without violating the group property, and without changing the code 
block length L, dimension of the code 77, and the code minimum distance d min . This means a more 
efficient code with the same parameters, i.e., 77 ,L and 7 rn in, would be available. The above discussions 
generalize to good codes with M > 2, including MDS codes. 

MIP entails that the coded symbols in each of A positions for A < 77 are independent of each other, 
each conveying log 2 (A7) bits of information about the coded message. This in turn means that each 
of the M constellation points occur with probability 1/M. Noting that the assignment of constellation 
points to M code alphabets is one-to-one and random, we can conclude that the constellation points 
in the A positions (selected by various code-words in a code satisfying MIP) are independent of each 
of other. 

Now let us assume the transmitted code-word is mapped to the constellation vector (sf, ..... sj 4 ), with 
each component selected from one of M possible Gaussian vectors with distribution J\f( 0, E s 1kxk)- 


'See j27| for definition and further details. 




The received vector is: 


( r i, i\l) = (sf,...,sf) + (z 1 ,...,z L ), (6) 

where (zi, zjf) is i.i.d. complex Gaussian distributed according to A/"(0, A^ 0 Ijr X ir)- Each component 
of the received vector (ri, ...,rt) is independently decoded to its nearest constellation point, resulting 
in a (hard decision) decoded vector denoted as (sf,...., s£). Now consider an error event £& with 
coordinates in any A < D positions. We have: 


£a = 


| 


B 


r h -si;|| > \\r h -s£||,...,||r, A -sfj\\ > ||r /A - 
z h \\ > \\z h + sf -s£||,..., ||z, A || > ||z /A + s£ -sf A 


)ll} 


j( s £ - s h) z h < - — 


sf — s 


B ||2 
h 
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I s *a s aII 


(7) 

( 8 ) 

(9) 


where z\ is the trasponse of vector z/. The components of the vector corresponding to such an error 
event, namely components of the vector (sf — sf,....,sf—sf), have Gaussian distribution (difference 
between two Gaussian vectors). We intend to show these components will tend to be independent of 
each other for large M. Indeed, as the underlying distributions are Gaussian, it suffices to show that 
the components of the difference vector, namely sf — sf and sf — sf for i ^ j are uncorrelated. 
Although: (1) constellation points are originally generated using an i.i.d. Gaussian distribution, and (2) 
it was shown that for A < D constellation points occurring in the A positions are independent of each 
other, the fact that different coordinates are populated with the same M realizations can be a source 
of dependency. Such a dependency occurs if two or more components of the vector (sf ,...., s f) take 
the same value. As the constellation points occur independently, each with probability 1/M, using the 
union bound, the probability of such repeated coordinates reduces at least with (A — 1)/M, which 
is negligible for large M. As the correlation coefficient between components of the difference vector 
is a weighted sum (over all occurrences of the difference vector) with such probability values as the 
weight factors, it follows that the contribution of such rare events to the correlation coefficient is 
negligible. It should be added that such partial contributions to correlation coefficient also depend on 
the square norm of the corresponding component of (sf, ....,s£) , which is finite and cannot change 
the above conclusion. In addition, M is independent of all other parameters contributing to such 
an asymptotical behaviour, which also confirms the validity of the earlier conclusion as M —>■ oo. 
Another scenario resulting in dependency occurs when a transmitted constellation point at a given 
coordinate, say l\, of (sf ,..... s/j is changed to sf in (sf ,..... s/^), and a different coordinate say 
1 2 in (sf,....,sf A ) also contains sf. The probability of this scenario can be expressed as F(Q 2 \ Gi) 
where, Q\ = {sf is transmitted and received in error as sf }, Q 2 = (sf = sf for some i / 1}. We 
have, 


P(ft Gi) = P(fe I stf, stf s,«) = P(& | sjj) = P( 62 ). 

This means the events Gi and G 2 are independent of each other. Again, using union bound, we can 
conclude that the probability of such event G 2 also reduces at least with (A — 1 )/M. Therefore, 
according to above arguments, events in equation [9] are independent of each other for large M: 
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( 10 ) 
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( 12 ) 


= E||v|| < <2 
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Vector v is complex Gaussian distributed according to A/”(0, 2 .E s Ik- x .k-). The probability density 
function for absolute value of such Gaussian distributed vector is: 

2 \ K-l 


P\M\(p) = 
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E S (K- 1)! \2E S 
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Defining Signal to Noise ratio, SNR = the above integral will be equal to (refer to [25]): 


P(£a) = I 2 


1 - 


SNR ( 1-SV 
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2 + SZVR ^ V k 

k =0 




(16) 


Asymptotic behaviour of P(£a) for large SNR can be expressed as: 

/ i / .'vv /-r \ ^ 

P(£a)~ 


1 - 


2 +a/vi? 


2 x SNR 


(17) 


For a FEC decoder with correction capability t = 1 j < D, the probability that the decoder be 

in error can be expressed as: 

p(s) = jr p(e A ) = p(s A ) + jr p(s A ). 

A=t+1 A=t+1 A=D+1 


( 18 ) 



























For A < Z), the arguments mentioned above concerning the independence of coded components will 
be valid, and we can write: 


D D 

E p < £ a>~ E 


A=t+1 A=t+1 
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On the other hand, as we are considering FECs with a large D, the term 
negligible and can be ignored. Overall, it follows that: 
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